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Abstract 

We consider a generahzed Brans-Dicke model in which the scalar field has a potential 
function and is also allowed to couple non-minimally with the matter sector. We assume a 
power law form for the potential and the coupling functions as the inputs of the model and 
show that acceleration of the universe can be realized for a constrained range of exponent 
of the potential function. We also argue that this accelerating phase is consistent with 
a large and positive Brans-Dicke parameter. In our analysis, the potential plays a more 
important role with respect to the coupling function in dynamics of the universe as the 
latter does not contribute to any of the relations characterizing evolution of scale factor 
of the universe and the scalar field. However, we will show that the coupling function 
is closely related to magnitude and direction of the energy transfer between matter and 
the scale field. We use this fact and some thermodynamic aspects of the model to put 
some constraints on the coupling function. In particular, we argue that the second law 
of thermodynamics constrains direction of the overall energy transfer. 

PACS Numbers: 04.50.Kd, 04.20. Cv, 95.36.+X 



1 Introduction 

Cosmological observations on expansion history of the universe can be interpreted as evidence 
either for existence of some exotic matter components or for modification of the gravitational 
theory. In the first route of interpretation one can take a mysterious cosmic fluid with suffi- 
ciently large and negative pressure, dubbed dark energy. In the second route, however, one 
attributes the accelerating expansion to a modification of general relativity. Such modified 
gravity models can be obtained in different ways. For instance, one can replace the Ricci 
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scalar in the Einstein- Hilbert action by some functions f{R) (for a review see, e.g., [1] and 
references therein), or by considering a scalar partner for the metric tensor for describing ge- 
ometry of spacetime, the so-called scalar-tensor gravity. The prototype of the latter approach 
is Brans-Dicke (BD) theory [2] which its original motivation was the search for a theory con- 
taining Mach's principle. As the simplest and best-studied generalization of general relativity, 
it is natural to think about the BD scalar field as a possible candidate for producing cosmic 
acceleration without invoking auxiliary fields or exotic matter systems. In fact, there have 
been many attempts to show that BD model can potentially explain the cosmic acceleration. 
It is shown that this theory can actually produce a non-decelerating expansion for low negative 
values of the BD parameter cu [3]. Unfortunately, this conflicts with the lower bound imposed 
on this parameter by solar system experiments [4] . Due to this difficulty, some authors propose 
modifications of the BD model such as introducing some potential functions for the scalar field 
[5], or considering a field-dependent BD parameter [6] without resolving the problem. 
In a general scalar-tensor theory there is a non-minimal coupling between the scalar field and 
Ricci scalar while the former minimally couples with the matter sector. In other terms, there 
is no an explicit coupling between the scalar field and matter systems in Jordan frame rep- 
resentation. In BD theory, in its original form, the motivation for such a minimal coupling 
was to keep the theory in accord with weak equivalence principle [2]. However, there has been 
recently a tendency in the literature to go a step further and consider a non-minimal coupling 
between the scalar field and matter systems as well by introducing an arbitrary function of the 
scalar field as a couphng function. In these models, the scalar field is regarded as a chameleon 
field as it can be heavy enough in the environment of the laboratory tests so that the local 
gravity constraints are satisfied. Meanwhile, it can be light enough in the low-density cosmo- 
logical environment to be considered as a candidate for dark energy. Such a chameleon-matter 
coupling was first introduced in BD model to achieve accelerating expansion of the universe for 
sufficiently large BD parameter [7] . It is shown that even though absolute value of uj is enlarged 
due to such a non-minimal coupling, its negative sign makes this theory remain inconsistent 
with observations. Later works apply non-minimal coupling to general scalar-tensor theories. 
For instance, [8] considers stability analysis and possibility of phantom crossing with an as- 
sumption that potential of the scalar field and the coupling function have power-law forms, or 
bouncing solutions and some cosmological tests are investigated in [9] . 

In the present work, our primary interest is to consider the possibility that such an anomalous 
coupling of matter systems can produce accelerating expansion of the universe for a sufficiently 
large positive BD parameter. As a first step, we will study conservation laws in such models 
in section 2. By writing the energy balance and the geodesic equations, we will show that 
how they are modified in terms of the introduced coupling function. It is also shown that 
the non-conservation of matter energy density means energy transfer between matter and the 
scalar field with a constant rate. In section 3, we study this model in a cosmological setting. 
We assume a power law form for the potential and the coupling function and show that ac- 
celerating expansion of the universe can be realized within a class of power law solutions of 
the field equations. In this class of solutions u can take positive and large values. We also 
discuss some thermodynamic aspects of the model and show that the overall energy transfer 
should be into the matter system if the second law of thermodynamics is to be fulfilled. We 
use this fact and some recent observations to tightly constrain the exponent of the coupling 
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function. 



2 The model 

We consider the action functional 

S^^J d'^x^icpR - ^5A^-V^0V,0 - 2V{<P) + 2/(0)L^} (1) 

where R is the Ricci scalar, is the BD scalar field, V{(p) and /(0) are some analytic functions. 
Here the matter Lagrangian density, denoted by L^, is coupled with via the function f{(f)). 
This allows a non-minimal coupling between the matter system and 0, thus the latter appears 
as a chameleon field. Taking f{(j)) — 1, we return to the BD action with a potential function 

Varying the action with respect to the metric g^i, and yields the field equations, given by, 

<f>G,. = + mTJZ (2) 

(2a; + 3)n<f) + 2{2V - V'<f)) = T"^/ - 2/'0L^ (3) 
where prime indicates differentiation with respect to 0, T"* = g^^TJJl and 

T^t = ^(V^</>V.0 - ^^^.V«0V» + (V^V.(/) - g,,aci>) - V{ci>)g^, (4) 

™ _ -2 (5(^/=ffL^) 

Due to the explicit coupling of the matter system with 0, the stress tensor is not divergence 
free. This can be seen by applying the Bianchi identities V^G^,/ = to (2), which leads to 

V^T- = {Lrng,. - T-) In / (6) 

As it is clear from (6), details of the energy exchange between matter and depends on the 
explicit form of the matter Lagrangian density L„i. Here we consider a perfect fluid energy- 
momentum tensor as a matter system with energy density pm and pressure Pm- 
There are different choices for the perfect fluid Lagrangian density which all of them leads to 
the same energy-momentum tensor and field equations in the context of general relativity [10] 
[11]. The two Lagrangian densities that have been widely used in the literature are Lm = Pm 
and Lm = —pm [12] [13] [14]. For a perfect fluid that does not couple explicitly to (i.e., for 
/(0) = 1), the two Lagrangian densities = p^ and = —pm are perfectly equivalent, 
as discussed in [13] [14]. However, in the model presented here the expression of enters 
explicitly the field equations and all results strongly depend on the choice of L^- In fact, it 
is shown that there is a strong debate about equivalency of different expressions attributed to 
the Lagrangian density of a coupled perfect fluid [15]. 
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Here we consider — Pm- In this case, although fluid elements follow geodesies of the 
background metric and there is no additional force, the matter is still non-conserved [16] 

Pm + 3H{^ + l)p^^-{^ + l)pmj (7) 

This has the solution 

Pm = Poa-'^^^'^ f-^^-"'^ (8) 

with po being an integration constant. This indicates that evolution of matter density strongly 
depends on the couphng function. When / = 1, (8) reduces to the standard evolution law for 
the matter energy density. However, as we have already stated, the coupling function in this 
framework does not affect a matter system with 7 = — 1 (a cosmological constant) as it is clear 
from (8). 



3 A cosmological setting 

Now we apply the above framework to a homogeneous and isotropic universe described by the 
metric (??). In a spatially flat case, the equations (2) and (3) become 

(2a; + 3)(0 + 3H^) - 2{2V - <I>V') = (1 - 3j)fpm + 2^<f>f'pm (11) 

To proceed further, we choose 1^(0) = V^cf)^^ and /(0) = /o0'^ in which li and I2 are constant 
parameters and Vq and /q are positive quantities. We assume power law forms for evolution of 
the scale factor and the scalar field 

a{t) = aot" (12) 

0(t) = (13) 

Inserting these ansatz functions together with (8) into (9)-(ll) gives the parameters n and m 
in terms of /i and I2. The latter parameters should be regarded as the inputs of the model 
described by (1) as they characterize the shape of the functions V{4>) and f{4>)- 
Since the left hand side of (9)-(ll) falls as t"'^, we arrive at the following relationships 

m{h - 1) = -2 (14) 

mih - 1) - (3n + mh) (7 + 1) = -2 (15) 
The equations (9) and (11), give then the consistency relations 

3n^ + 3mn - \m'uj - H</>o^"' - /o">oao '^^^'Vo ^''''^'^ - (16) 
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(2a; + 3)(m(m - 1) + 3mn) - 2(2 - h)Vo(j)'o'~^ + [(3 - 2^2)7 - IJ/tTVoao "^"'^'Vo ^""'""'^ = (17) 
Combining (14) and (15) gives 

" = 3(,+ l)ft-l) (^**) 

For 7 7^ 0, the conditions for the universe to be expanding (n > 0) and accelerating {n > 1) 
are the following^ 

n > l2> -/1/7 Zi > 1 

k < -hh h < 1 (19) 

n>l ^ k>hh-l) + ^{h-3) k>l 
z Z7 

i2<^(ii-l) + 7^(/i-3) li<l (20) 

2 27 

As an illustration, let us consider the potential V{(f)) = Vo0^ in a radiation-dominated universe. 
In this case, 7 = 1/3 and expansion of the universe implies I2 > —6. On the other hand, the 
success of big bang nucleosynthesis gives us a strong evidence of the radiation-dominated 
decelerated phase which leads to Z2 < 0. Thus a scalar field with a quadratic potential which 
couples with radiation requires that — 6 < ^2 < 0. Note that the relation (14) indicates that for 
/i > 1 the scalar field is a decreasing function of time (m < 0). This implies that the coupling 
function / oc t"*'^ increases with time or the matter-chameleon coupling gets stronger as the 
universe expands. 

This feature is changed in a matter-dominated universe. For 7 = 0, the expression (18) reduces 
to 

This relation indicates that n > 1 can not be reahzed if h < 0. On the other hand, when 
Zi > accelerating expansion of the universe requires that 

l<h<3 (22) 

The parameter I2 does not enter the above condition and it seems that the scalar field effectively 
decouples from the matter system in a matter-dominated universe. Exploring the relations (14) 
and (15) reveals that I2 also disappears from these relation when we set 7 = 0. The same is 
true for the consistency relations (16) and (17). 

As it is evident, the above conditions (19), (20) and (22) do not put any constraint on the BD 
parameter ou. To obtain such a constraint, one should consider the right hand side of (10) as 
a function 



'''We assume that the matter system satisfies the weak energy condition 7 + 1 > 0. 
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and establish the condition for having C > 0. For 7 = 0, the condition is 

At this point, attention is focused on the question that whether it is possible to have accelerated 
expansion for a large positive uj. To answer the question, let us consider the above expression 
for u » 1, 

200 ti - 1 

If the left hand side effectively takes a negative sign then the inequality can not be satisfied 
for a positive BD parameter. To make an estimation, we first note that (p^-^ ~ G ~ and 
Po^o^ ~ M with G, Mp and M being respectively the gravitational constant, the Planck mass 
and the total mass content in the universe. Then we can rewrite (25) in the form 

It is an observational fact that matter density of the universe is of the same order of the critical 
density pc = SH'^M^ /8n with H being the Hubble parameter. This leads to a relation between 
M and H such that M/M^ ~ R, which R ~ H'^ is the Hubble radius^. Substituting this 
result into (26), gives 

M?'-»(V„ -^^)-^^>^^^ (27) 

Using R ~ lO^^(meter) and M ~ 10^^(meter)~^^, one can see that the term containing Vq is a 
leading term on the left hand side of (27) unless li takes values near the lower bound in (22), 
namely /i ~ 1. In other words, when li is not close to unity one can write cj < M^('i~^^Vo/m^ 
which indicates that the BD parameter can take positive large values. Note that there is no 
need for fine-tuning of the constant Vq since the upper bound given by the latter relation 
is sufficiently large due to the appearance of Mp. It can be easily checked that for /i 1, 
R^i IM^^~^ ^ R 10^^{meter) and the left hand side of (27) take a negative sign. In this 
case, the inequality can not be satisfied for u > 0. It is interesting to note that the coupling 
function, or the parameter I2, does not play any role in the above argument and it is only the 
potential function that is relevant. As we have already seen, this is also the case (for 7 = 0) 
when one enters the arguments concerning the expansion of the universe in matter-dominated 
era. 

Despite the irrelevant role of the coupling function in the latter arguments, it affects the 
dynamics of the matter energy density. Combining (12) and (13) with (8) gives 

Pm oc a-3(^+^)+^ (28) 

where 

771 

£ = ^2(7+1) (29) 

n 



■I- One can say that the radius of the universe coincides with its Schwarzschild radius 2GM [20]. 
^We use the unit system in which U = c = 1. 
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In fact, the parameter e measures the modification of matter expansion rate due to its in- 
teraction with the scalar field. The relation (28) states that when e > matter is created 
and energy is constantly injecting into the matter so that the latter will dilute more slowly 
compared to its standard evolution oc a^'^^'~'^^K Similarly, when £ < the reverse is true, 
namely that matter is annihilated and direction of the energy transfer is outside of the matter 
system so that the rate of the dilution is faster than the standard one. 
We may use (15) to write (29) in the form 

n - ^^^(^ + 1) (30) 

3/2(7 + 1)^-^(7^2 + 1) ^ ^ 

When m = 0, (j) takes a constant configuration which is given by 0o ~ and the model (1) 
reduces to Einstein gravity with a cosmological term. This case corresponds to £ = and the 
solution (30) reduces to n = 2/3(7 + 1), the case in the standard cosmology. 
For 7 = 0, the expression (30) takes the form 

^ - TT^ (31) 
h 

One can use this relation to constrain the parameters e and I2. We first note that from the 
modified matter expansion rate (28) we should expect \e\ « 1 since so far there has been no 
report from observations about an anomalous expansion rate. In other terms, if there exits 
any deviation from the standard evolution of matter density it must be very small. Taking e 
as an independent parameter, one can compare the modified expansion rate with observations 
and perform data fitting to estimate \e\. In fact, it is shown that cosmological observations are 
consistent with \e\ < 0.1 [17] [18] [19]. This result puts an upper bound on absolute value of 
the parameter I2. The condition for accelerating expansion inferred by (31) is e/l2 > 1 which 
leads to I/2I < 0.1. 

On the other hand, the expression (31) implies that the parameters e and I2 should have the 
same sign in order that n > 1. Thus direction of the energy transfer between the matter 
system and the scalar field is characterized either by e or /2- Following [18], we argue that 
£ > as required by the second law of thermodynamics. To do this, we should investigate some 
thermodynamic features of the matter-chameleon coupling described by (1). A thermodynamic 
description of a perfect fluid matter system requires the knowledge of the particle flux A^" = 
nu°' and the entropy flux S°' = nau'^ where n = N/a^ and a = S/N are, respectively, the 
concentration and the specific entropy (per particle) of the created or annihilated particles. 
Since the energy density of matter is given by pm = nM with M being the mass of each 
particle, the appearance of the extra term in the energy balance equation (7) means that this 
extra-change of Pm can be attributed to a change of n or M. Here we assume that the mass 
of each matter particle remains constant and the extra term in the energy balance equation 
only leads to a change of the number density n. In this case, the equations (7) can be written 
as (7 = 0) 

h + 3Hn = nV (32) 

where 

r.-i = .^ (33) 
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is the decay rate. We also assume that the overall energy transfer is an adiabatic processes in 
which matter particles are continuously created or annihilated while the specific entropy per 
particle remains constant during the whole processes {& — 0) [21]. This means that 



S _N 



(34) 



From n oc a^^+^, we can sec that the total number of particles scales as oc a^, and 
the second law of thermodynamics 5" > imphes that £ > 0. This conclusion can also 
be drawn by (33) since F > requires that £ > in an expanding universe. Hence, the 
chameleon scalar field should suffer energy reduction and the matter system should gain 
energy during expansion of the universe if the second law of thermodynamics is to be ful- 
filled. 



4 Conclusions 

In this work we have studied some features of the generalized BD model in which the scalar 
field is allowed to couple non-minimally with matter sector. The matter expansion law and the 
geodesic equation have modified due to this non-minimal coupling. The modification depends 
on the choice of L^, Lagrangian density of the matter system. For instance, in our choice 
matter conservation equation is modified while there is no extra force in the geodesic equation 
since = leads to vanishing of the first term on the right hand side of the equation (??). 
We have found power law solutions for the scale factor and the scalar field. In this class of 
solutions, when 7 = accelerating expansion of the universe can be realized for 1 < li < 3. 
Our analysis also indicates that this accelerating phase is consistent with a positive and large 
BD parameter. 

The fact that the parameter I2 does not contribute to the condition for cosmic speed-up in the 
case of 7 = 0, may be regarded as a signal for an irrelevant role of the function /(0) in the 
matter-dominated era. However, it is evident from (8) that /(</>) is important in the evolution of 
matter density. We have reformulated the matter expansion law as p„j oc a~^^'^'^^^'^^ where the 
parameter e characterizes both magnitude and direction of the energy transfer. We have argued 
that such a reformulation has two important consequences. First, since the magnitude of e 
can be fixed by observation (as it is recently suggested that \e\ < 0.1) accelerating expansion 
which requires that I2 < e sets an upper bound on the absolute value of 12- Second, the sign 
of e is restricted by the second law of thermodynamics to assume only positive values so that 
the direction of energy transfer is into the matter system. This constrains the exponent of the 
coupling function to take values within the range < Z2 < 0.1. 
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